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DESIGN OF AN OIL SQUEEZE FILINA DAMPER BEARING FOR A MULTIMASS 
FLEXIBLE-ROTOR BEARING SYSTEM 

by Robert E. Cunningham, Edgar J. Gunter, Jr. /and David P. Fleming 

Lewis Research Center 

SUMMARY 

Results of a single-mass flexible-rotor analysis were used to optimize the stiffness 
and damping of a flexible, damped support for a symmetric five -mass rotor operating 
through and above its first bending critical speed. The flexible support attenuates the 
rotor amplitudes and forces transmitted to the ball bearings. 

The stiffness and damping values obtained from the single-mass analysis were then 
used in the design of an oil squeeze-film damper and its centering spring. Short bearing 
lubrication theory was used. A cavitated oil film was assumed to exist in the nonrotating 
damper journal. 

The damper design was verified by an unbalance response computer program. At 
the first critical speed, rotor amplitudes were reduced by a factor of 16 and bearing 
loads reduced by a factor of 36 compared with the same rotor with rigidly supported 
bearings. Amplitudes and forces at higher speeds were also reduced substantially. 


INTRODUCTION 

In many rotor design applications such as in turbojet engines, compressors, or tur- 
bines, the rotor experiences high vibrational amplitudes resulting in large forces being 
transmitted to the bearings and the support structure. These high vibrational responses 
may be due to several causes and may be roughly grouped under the headings synchro- 
nous and nonsynchronous response. The forces and amplitudes for synchronous response 
are usually associated with unbalance forces in the rotor. This unbalance may be a re- 
sult of either the manufacturing process and/or the assembly of the components. Even 
if a rotor is well balanced initially, the balance degrades with use. Thermal gradients 

* Professor of Mechanical Engineering, University of Virginia, Charlottesville, 
Virginia. 


can cause the shaft to warp. Erosion of compressor and/or turbine blades can alter 
the balance of the rotor assembly. Therefore, in the design of a turbomachine provis- 
ions should be made so that the increase of unbalance with engine operation will not load 
the bearings excessively or cause large rotor amplitudes. 

Another serious problem related to high-speed turbomachinery is the occurrence 
of nonsynchronous, self-excited whirl motion. This is commonly associated with fluid 
film bearings, but can also be caused by rotor internal friction (ref. 1) or variable 
aerodynamic loading (ref. 2). 

Theoretical studies conducted by one of the authors (ref. 3) indicated that problems 
of both self -excited rotor instability and high vibrational (synchronous) response can 
often be greatly alleviated by a properly designed damping system at the rotor supports. 
Flexible damped rotor supports may be used to 

(1) Reduce the forces transmitted through the bearings and foundation 

(2) Reduce the amplitudes of motion of the rotor that could result in rubbing and 

excessive wear of close fitting components 

(3) Permit smooth operation through critical speeds 

(4) Protect the machine from sudden buildup of unbalance forces due to compressor 

or turbine blade loss 

(5) Protect the machine from potentially destructive, self-excited instability 
This report will consider only steady, synchronous motion; that is, effects (1) to (3). 

Damping may be achieved by various mechanisms, such as Coulomb friction, vis- 
coelastic materials, and viscous dampers, which may utilize either compressible or 
incompressible fluids. This investigation will be concerned with the damping character- 
istics obtained for the incompressible-oil squeeze-film damper. Such a damper appears 
in figure 1, which shows a rolling-element bearing mounted in an oil squeeze-film 
damper. The annulus between the outside diameter of the ball bearing housing and the 
damper housing inside diameter is filled with oil. The orbital motion or precession of 
the ball bearing housing in the damping fluid generates a hydrodynamic pressure. This 
particular type of damper is now being used in several production aircraft turbojet en- 
gines and in other types of high speed turbomachinery. 

This investigation was conducted to (1) examine the influence of flexible damped 
supports on rotor amplitudes and forces transmitted over a specific operating speed 
range, (2) show how single-mass rotor theory can be used to design a support system 
for a multimass rotor operating below the second bending critical speed, and (3) demon- 
strate design procedure for an oil squeeze-film damper. 


SYMBOLS 

A amplification factor at rotor critical speed, dimensionless 
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B b 

! d 

B d 

B s 

b i 

B 2 

C 

C* 


c 


r 


E 

e 

e. 


At 


K>* t 


H 

h 

K 

K c 

K d 

K d 

K s 


support amplitude, m (in. ) 

rotor amplitude, m (in. ) 

amplitude at the critical speed, m (in. ) 

damping ratio, dimensionless 

bearing damping, N sec/m (lb sec/in. ) 

squeeze film damping, N sec/m (lb sec/in.) 

squeeze film damping coefficient, dimensionless 

shaft damping, N sec/m (lb sec/in. ) 

support damping, N sec/m (lb sec/in. ) 

effective rotor-bearing system damping, N sec/m (lb sec/in.) 
basic dynamic capacity of ball bearing, N (lb) 
adjusted dynamic capacity of ball bearing, N (lb) 

one -half total clearance in damper bearing, also radial damper clearance, 
m (in. ) 

shaft diameter, m (in.) 

O 9 

elastic modulus, N/m (lb/in. ) 
bearing eccentricity, m (in.) 
rotor mass eccentricity, m (in.) 
bearing force, N (lb) 

radial and tangential bearing forces, N (lb) 
force transmitted to foundation, N (lb) 
gravitational constant, m/sec (in. /sec ) 
ball bearing life, hr 
fluid film thickness, m (in. ) 
stiffness ratio = Kj/Kg, dimensionless 
bearing stiffness, N/m (lb/in.) 
centering spring stiffness, N/m (lb/in.) 
squeeze-film stiffness, N/m (lb/in. ) 
squeeze-film stiffness, dimensionless 
shaft stiffness, N/m (lb/in. ) 
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Kjl support stiffness, N/m (Ib/in. ) 

K 2 effective rotor -bearing system stiffness, N/m (lb/in. ) 
L bearing or shaft length, m (in. ) 

M mass ratio = M^/M 2 

Mj total support or bearing housing mass, kg (lb sec /in.) 
M 2 rotor mass, kg (lb sec / in.) 

N rotor speed, rpm 

9 

p pressure, N/m (psi) 

R damper bearing radius, m (in. ) 

t time, sec 

U unbalance, kg m (oz in. ) 

u circumferential journal velocity, m/ sec (in. / sec) 

z axial coordinate, m (in. ) 

e eccentricity ratio = e/c , dimensionless 

Q, 9 ' angular displacements, rad 

/x absolute viscosity, N sec/mr (lb sec/in. ) 

journal precession rate, rad/sec 

u> angular velocity, rad/sec 

o> rigid support critical speed, rad/ sec 

c r 


INFLUENCE OF DAMPER SUPPORT ON SINGLE- MASS FLEXIBLE 
ROTOR ON ROLLING-ELEMENT BEARINGS 

Successful operation of a high-speed rotor requires careful balancing to minimize 
the amplitude of motion and the forces transmitted through the bearings. Obviously, 
the magnitude of the forces transmitted will have a significant influence on the life of a 
rolling-element bearing. 


Example 1 - Load, Life Relation for Rolling -Element Bearing 

Consider a flexible rotor operating at a speed of 30 000 rpm and supported on rigidly 
mounted, 204 size, 20 -millimeter bore, ballbearings. The basic load rating (see 
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ref. 4) of the bearing C is 9800 newtons (2210 lb). If the bearing is to have a reliabil- 
ity of 99 percent, the load reduction factor as given by Harris (ref. 5) is C'/C = 0.61. 
For a design life of 2000 hours the permissible bearing load is 



(or 344 lb) 


If the rotor has an unbalance of U acting at the bearing, the unbalance force trans- 
mitted is given by 


F, = U(oj 2 ) = U — (1) 
D 30 

for F = 1530 newtons (344 lb), the unbalance U is 15. 5 gram-centimeters (0. 216 
oz in.). Note that, because this is a flexible rotor, the effect of static unbalance may 
be magnified by rotor bending. Thus, the unbalance may need to be considerably lower 
than calculated to keep the bearing load to a permissible level. Figure 2, calculated 
from the equations in reference 3, shows the forces acting on the bearings over a range 
of speed. Whenever the rotor speed is less than 1. 4 times the critical speed, the bear- 
ing force is greater than predicted by the simple equation (eq. (1)). At the critical 
speed, for this lightly damped rotor, the bearing force is 10 times that for a rigid rotor. 
If the unbalance is doubled, 31. 0 gram -centimeters (0. 432 oz in.), the calculated ro- 
tating load increases F to 3060 newtons (688 lb). The basic load rating equation shows 
that doubling the bearing load will cause a reduction of bearing life by a factor of 8. 

Thus 


H = 250 hr 

From the preceding example, it can be seen that with a rigid support the unbalance can 
cause large forces to be transmitted through the bearing. If, however, the rotor is 
mounted on flexible, damped supports, then the forces transmitted through the bearings 
may be attenuated by the proper selection of the support damping value. 

Figure 3 shows schematically a single-mass rotor mounted in elastic damped sup- 
ports. Figure 4 (from ref. 3) represents the rotor amplitude against rotor speed for 
various values of support damping for this system. In figure 4 rotor speed has been 
normalized with respect to the rigid support critical speed. If the bearing supports have 
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no damping, the rotor then will have a very high response at the critical speed, which 
is 3/10 of the rigid support critical speed. 

As the damping is increased in the flexible support, the rotor amplitude at the crit- 
ical speed is diminished until, at the optimum value of dimensionless support damping 
(B = 5. 0), the peak amplitude at the critical speed completely disappears. However, if 
the support damping value is increased beyond the optimum value, the amplitude of mo- 
tion will increase at the rigid support critical speed. For example, a damping value 
of 50 is excessive and has caused the support to lock up. Thus the behavior approaches 
that of a rigidly supported rotor. Figure 4 has been plotted for an amplification factor 
A of 10. The amplification factor is the ratio of the rotor amplitude at the rotor crit- 
ical speed to the rotor unbalance eccentricity, that is, A = a cr /e^. A value of 10 rep- 
resents moderately light damping. Reference 3 points out that optimum support damp- 
ing is virtually independent of A for A > 10; thus the information of reference 3 is 
applicable to a wide range of rotors. 

The amplification factor (as shown in ref. 3) can also be expressed as 


w cr®2 


where Bg is the equivalent rotor and bearing damping given by 


B 2 = 


K s B b 


+ K / + ("V 2 


+ B, 


The critical speed o> cr is calculated from 



(3) 


(4) 


where Kg is the equivalent rotor and bearing stiffness given by 

K W K s * Kfr> + K^V 2 

2 (Kb + Kg) 2 + (wiy 2 

When the bearing damping is small (as for rolling -element bearings), the equivalent 
stiffness and damping reduce to 
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( 6 ) 


, «b K s 
Kb + K s 


and 


b 2 = 




+ B„ 


<®b + V 


( 7 ) 


The preceding equations (2) to (7) are given in reference 3. 

The data in figure 4, which are for a low mass ratio system (M = 0. 1), can be ap- 
plied to the single-mass rotor shown in figure 5. Example 2 will illustrate how to de- 
termine the support stiffness and damping required to prevent the excitation of the rotor 
first bending critical speed. 


Example 2 - Magnitude of Damping Required to Attenuate Rotor Amplitude at 
First Critical Speed for Low Support Stiffness and Mass Ratio 

The rotor of figure 5 is modeled as a mass of 2. 42 kilograms (0. 0138 lb-sec 2 /in. ) 
on a massless elastic shaft. The spring rate for this rotor is calculated from 

3ffED^ 

V _ s / 


where 


E =2. 1X10* 1 N/m 2 


(or 30x10® psi) 


D g = 0.0254 m (or 1 in. ) 
L = 0. 48 m (or 19 in.) 


K = 18. 2X10® N/m 


(or 10 300 lb/in. ) 


Thus 


The rotor is mounted in ball bearings which are assumed to have a much higher stiff- 
ness than that calculated for the shaft. Thus the critical speed on rigid supports can be 
calculated from 



oj = 8280 rpm (867 rad/sec) 

The amplification factor A for this rotor on rigid supports is assumed to be 10. 
Solving equation (2) for Bg results in 


B„ = = 210 N‘ sec/m (or 1. 20 lb - sec/in. ) 

cr 


The optimum damping ratio B for a minimum response over the speed range was 5 
(see fig. 4). Solving for the required support damping results in 

B^ = B x Bg = 1050 N sec/m (or 6. 00 lb sec/in. ) 

The corresponding support stiffness for a stiffness ratio K of 0. 1 is ^ 

K t = K x K 2 = K x K g = 0. 1 K s = 1. 82X10 5 N/m (or 1030 lb/in. ) ; 

l i 

In an actual design a support with such a soft spring rate would be susceptible to shock 
and self -excited whirl instabilities (ref. 3). A support having a higher stiffness may be 
desirable; therefore, a different value of damping would be necessary to minimize the 
rotor amplitudes and forces transmitted. 

A further consideration in the choice of support stiffness is the mass ratio M, 
which is the ratio of bearing housing mass Mj to rotor mass Mg. Figure 6 shows 
that the mass ratio M, if one has a choice of stiffness ratios, should be as low as pos- 
sible to minimize the rotor amplitude. For a given mass ratio, a tuned system produces 
near -minimum amplitude. A tuned system is a support system for which the stiffness 
ratio K equals the mass ratio M. Since for some rotor -bearing systems it is difficult 
to make the bearing and its housing much lighter than the rotor, the following example 
will use mass and stiffness ratios of unity, that is, M = K = 1. 
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Example 3 - Damping Required to Minimize Rotor Amplitude When K = M = 1 

Figure 7 shows optimum damping ratios for a tuned system and the resultant mini- 
mum rotor amplitudes. For a value of M = K = 1, the optimum damping ratio B is 14. 
Rotor amplitudes over a speed range for several damping ratios are shown in figure 8. 
The rotor -bearing damping B 2 was previously calculated as 210 newton- seconds per 
meter (1. 20 lb-sec/m). Thus the support damping is given by 

= B X B 2 = 2940 N* sec/m (or 16. 8 lb ♦ sec/in. ) 

From figures 7 and 8 for K = M = 1 and the optimum damping ratio of 14, the 
maximum rotor amplitude will be 1.6 times the unbalance eccentricity. This represents 
a reduction in amplitude to one -sixth that of a rotor running on rigid supports that has 
an amplification factor of 10. 

In the rigidly mounted bearing model considered in the sample ball bearing calcula- 
tion (example 1), the ratio of the transmitted force to the rotating unbalance force ranged 
up to 10. In the proper design of a flexibly mounted rotor with damping, the dynamic 
transmissibility and the bearing forces should be considerably less. With insufficient 
damping the dynamic transmissibility may exceed the rigid support value at the rotor 
critical speed (see fig. 9 for B = 0. 01). In this instance the force transmitted through 
the bearings also exceeds the rigid support value (fig. 10). The problem to be consid- 
ered now will be the selection of a value of damping to use in the support to minimize 
bearing forces and forces transmitted to the foundation. Figure 9 is a plot of the dy- 
namic transmissibility over a speed range for various values of damping. The damping 
value chosen to minimize /rotor amplitude (B = 14) is near -optimum for minimizing 
transmissibility, although B = 10 results in a slightly lower maximum foundation force. 
Figure 10 shows that B » 14 minimizes the force transmitted through the bearings. 

When the rotor is operating above 1. 4 times the rigid support critical speed, support 
damping decreases the forces transmitted through the bearings, although increasing the 
damping at high speeds (>1. 8 u> cr ) will cause increased loads to be transmitted to the 
foundation. The proper design of the damper then must take into consideration the 
forces and amplitudes of motion to be permitted throughout the operating speed range. 

In addition to these steady-state conditions, the proper damper design must also take 
into consideration the damping necessary to insure stability and transient operation due 
to shock and other suddenly applied loads (such as those that might occur due to the loss 
of a blade in the compressor or turbine). For the rotor considered here a damping ra- 
tio B of 10 to 14 should be satisfactory. Since the rotor -bearing damping B 2 is 
210 newton -seconds per meter (1. 20 lb sec/in. ), the required support damping is 10 to 
14 times this, or 2100 to 2940 newton -seconds per meter (12. 0 to 16. 8 lb-sec/in. ). 
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DAMPER SUPPORT REQUIRED FOR A MULTEMASS ROTOR OPERATING 
BELOW THE SECOND BENDING CRITICAL SPEED 

In contrast to the singlemass rotor, a multimass rotor has more than one bending 
critical speed. Thus, in general, singlemass rotor data are not adequate for the design 
of a damper support for a multimass rotor. However, if the maximum service speed 
is below the second bending critical speed, singlemass theory may be entirely adequate, 
as will now be shown. 

The rotor shown in figure 11 is designed to simulate the flexible rotor of a small, 
lightweight turbocompressor which is required to operate above its first bending critical 
speed and below the second bending critical speed. Oil squeeze film dampers are to be 
used at the bearing supports to attenuate rotor motion. Two single -row, deep -groove 
ball bearings, series 204, will be used to support the rotor. Each bearing support 
housings has a mass of 1. 21 kilograms (0. 0069 lb sec/in. 2 ). A cantilevered centering 
spring supports the ball bearing housing; the spring rate can be chosen to complement 
that of the squeeze film. Oil is supplied to the damper from a circumferential groove. 
Two piston rings located in circumferential grooves along with metering orifices control 
the flow of damping oil from the bearing ends. These features are shown schematically 
in figure 1. 

Critical speeds were calculated for this rotor on rigidly supported bearings by the 
critical speed computer program of reference 6. A bearing stiffness of 65. 5 meganew- 
tons per square meter (375 000 lb/in. ) was assumed. The first three critical speeds 
and associated mode shapes are shown in figure 12. The first bending critical speed 
was calculated as 8280 rpm (867 rad/sec). 

The shaft and bearing span of this five-mass rotor are identical to those of the 
single-mass rotor of example 2. Thus, the five-mass rotor will have the same stiffness 
as the single-mass rotor. For calculation purposes, the five masses may be replaced 
with a single mass at the rotor center, which results in the same critical speed. This 
equivalent single mass may be calculated from 

m 2 = — = — MNy/m2 = 2. 42 kg (or 0. 0138 lb ■ sec 2 /in. ) 
o> 2 r (867 rad/sec) 2 

The actual rotor mass is 5. 68 kilograms (0. 0325 lb sec 2 /in. ). 

The equivalent mass is identical to that of the single mass rotor. The ratio of sup- 
port mass to rotor mass is 
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M = — = . 1 

M 2 2. 42 

The supports will be designed for a tuned system, that is, K * M. This being the case, 
the rotor properties, support stiffness, and mass ratios are the same as those for the 
single-mass rotor, and all of the previously calculated values may be used without 
change. Thus the required damping ratio B is 10 to 14, and the support damping 
is 2100 to 2940 newton -seconds per meter (12. 0 to 16. 8 lb sec/in. ) or one-half of this 
for each of the two supports. 


CHARACTERISTICS OF THE SQUEEZE-FILM DAMPER 

A squeeze -film damper is shown schematically in figure 13. It consists of a cylin- 
drical journal, which is prevented from rotating in a cylindrical bearing. The journal 
center is assumed to make a circular orbit about the bearing center. 

The Reynolds lubrication equation for incompressible flow is given in reference 7 as 

3e\ + i-/i>!8E\ = 6 3ihu2 +12 3h (W 

r 2 00 00/ 0z yju 0z J R 00 0t 


If the bearing is considered to be very short relative to its diameter, the pressure- 
induced flow (h 3 / p)(0 p/00) in the circumferential direction is negligible compared with 
the shear flow in the circumferential direction. It is convenient at this point to intro- 
duce rotating coordinates according to 


0’ = 0 - cot 


( 10 ) 


where u> is the rate of precession. 

The Reynolds equation now becomes 



Assuming that the bearing and journal axes are parallel, the film thickness h is given 

by 


h = c r (l + e cos 0') 


( 12 ) 
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where the eccentricity ratio e is defined by 


e 


e 



(13) 


Since h does not depend on z, the Reynolds equation may be integrated directly with 
the boundary conditions P(0',O) = P(0 f , L) = 0 to give 

P(0',z) = (Lz _ z 2) , 

h 3 30' 


It is common practice in lubrication analyses to neglect negative pressures since* in 
general, a fluid cannot sustain a tensile stress. If it is assumed that cavitation will 
occur, the preceding expression may be integrated over the area of positive pressure to 
give a radial restoring force and a tangential force that opposes the journal 
motion. The results, which are derived in detail in reference 8, are 


F 


r 


, 2jiRL 3 q?€ 2 



3 

nfiRh we 

2c 2 r (l -Sf 2 


(15) 


(16) 


These results can also be given in terms of stiffness and damping coefficients: 




2uRL 3 €0) 



(17) 


and 


B d 



(18) 
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These coefficients may be expressed in dimensionless form as 


and 



If n 3 

K d c r 


l±w RL' 



( 19 ) 



BjC 3 

dr 

/i.L 3 R 



( 20 ) 


Both and B d are functions only of the eccentricity ratio and are plotted as 
such in figure 14. 


DESIGN OF SQUEEZE-FILM DAMPER BEARING FOR MULTIMASS FLEXIBLE 
ROTOR TO OPERATE THROUGH THE FIRST BENDING CRITICAL SPEED 

The magnitude of damping required has been determined, and it is now necessary 
to design the damper bearing to produce this amount of damping. The expression for 
B d ( e q- (18)) shows that damping is a function of damper clearance, length, radius, ec- 
centricity ratio, and oil viscosity. Generally, the diameter of the damper housing is 
dictated by the ball bearing outside diameter. The same oil is usually used for both the 
damper and the ball bearing; thus, the viscosity of damper oil is fixed. It remains, 
however, to select values of the radial clearance c r and length L. 

Since film pressure and consequently film stiffness increase rapidly with increasing 
eccentricity ratios, it is not desirable to operate at very large eccentricity ratios be- 
cause the film stiffness will make the overall support stiffness too large. A maximum 
eccentricity ratio of e = 0. 4 at the first critical speed was chosen for this damper de- 
sign. For that value a dimensionless damping B d of 2. 04 is obtained from figure 14. 

The next parameter to determine is the damper clearance c r> Figure 15 (from 
ref. 3) shows that, for an optimally damped system (B = 10 to 20), the maximum support 
amplitude is about equal to the mass eccentricity of the rotor, that is, the displacement 
of the rotor center of gravity from the geometrical center. The five-mass rotor of fig- 
ure 11 is expected to have a maximum unbalance of 29 gram-centimeters (0. 4 oz in. ) 
distributed fairly uniformly over the five masses. This corresponds to a mass eccen- 
tricity of 0. 05 millimeter (0. 002 in. ). If the maximum damper eccentricity ratio is 
0. 4, the damper clearance must be 0. 05/0. 4 = 0. 13 millimeter (0. 002/0. 04 = 0.005 in. ). 
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It remains only to determine the damper length. The design of figure 1, with a 
circumferential oil supply groove divides the damper into two equal halves. Each half 
functions separately and thus provides one-fourth of the total required support damping. 
The total support damping needed was previously determined to be 2100 to 2940 newton 
seconds per meter (12. 0 to 16. 8 lb sec/in. ). The rotor is expected to operate with con- 
siderably less than the maximum unbalance most of the time; the maximum of 29 gram- 
centimeters (0. 4 oz-in. ) represents a degraded value after considerable service time. 
Less unbalance means lower amplitudes and, thus, lower stiffness of the squeeze film 
(eq. (17)). To maintain squeeze-film stiffness as much as possible at lower eccentric- 
ities, the damper will be sized for damping near the top of the range for an eccentricity 
ratio of 0. 4. Thus, the design value will be 2800 newton-seconds per meter 
(16 lb sec/in. ). Damper length can now be determined from 

L = c 


where 



c r = 0. 13 mm (or 0. 005 in. ) 

^140° F = °- 0119 N sec / m2 (° r 1- 73X10 -6 lb sec/in. 2 ) 

R = 39. 6 mm (or 1. 56 in. ) 

B^ = 2800 N-sec/m (or 16 lb-sec/in.) 

B d = 2.04 

L = 11. 4 mm (or 0. 45 in. ) 

The net radial stiffness of the cavitated squeeze film can now be calculated. For 
e = 0. 4 a value of dimensionless stiffness K d = 1. 13 is obtained from figure 14. A 
value of film stiffness can now be calculated from 


K d = 


K d p(oRL 3 


c 


3 

r 
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where co is the damper angular precession speed. At the first critical speed of 
867 rad/ sec (8280 rpm), = 338 kilonewtons per meter (1930 lb/in. ). The total 
damper spring rate is four times this value or 1350 kilonewtons per meter (7720 lb/in. ). 
As was stated before, the overall support stiffness is the sum of the centering spring 
stiffness and the damper stiffness since the two springs act in parallel. The centering 
spring stiffness may now be chosen to provide optimum stiffness at the first critical 
speed. For K = 1 total support stiffness = 1820 kilonewton-meters (10 300 lb/in.). 
Thus, the stiffness K c of each centering spring must be 1/2(1820 - 1350) = 235 
kilonewton-meters (1290 lb/in. ). This is a very soft spring and will need to be pre- 
loaded to center the rotor at low speeds. 

To determine how effective this damped flexible support is in attenuating the rotor 
amplitude and bearing forces, an unbalance response computer program was used to 
produce figures 16 to 18. The program is that of reference 9 which treats nonaxisym- 
metric rotor supports and nonlinear stiffness and damping in the squeeze film. 

Figure 16 shows rotor amplitude at midspan, and figure 17 shows forces that would 
be transmitted to the ball bearing for both a rigid and a flexible damped support. The 
long-dash curves of these figures represent the rotor operating with good balance; the 
total imbalance is 7 gram-centimeters (0. 1 oz in. ). The corresponding mass eccentric- 
ity is 0. 013 millimeter (0. 0005 in. ). The first critical speed has shifted upward to 
about 9000 rpm. Figure 16 shows that the rotor amplitude is about three times the 
mass eccentricity at the first critical speed. The amplitude then drops with increasing 
speed and remains low out to 30 000 rpm. Amplitude then increases with speed, reach- 
ing a maximum of four times the mass eccentricity at 39 000 rpm. Figure 17 shows 
that the bearing experiences virtually no force buildup due to the critical speed. At 
9000 rpm the bearing force is only 14 newtons (3 lb). This is only one-half the force 
that would be experienced by a bearing on a rigidly supported rigid rotor (for which 
F^ = (l/2)M2e^w^). Bearing force generally rises with speed to 39 000 rpm and then 
drops off. 

Now consider the case where the initial balance at assembly has degraded for any 
one or a combination of the reasons mentioned in the INTRODUCTION. Instead of a 
total unbalance of 7 gram -centimeters (0. 1 oz-in. ), let us assume the unbalance now 
is 29 gram -centimeters (0. 4 oz-in.) resulting in a mass eccentricity of 0. 05 millime- 
ter (0.002 in. ). The solid curves of figures 16 and 17 show results for the rotor on 
squeeze -film supports, and the short-dash curves show results for the rotor on rigid 
supports. At the first critical speed the center amplitude of the flexibly supported rotor 
is about three times the mass eccentricity, as with the lesser unbalance. Over the en- 
tire speed range the amplitude changes nearly the same as for the low unbalance. How- 
ever, the peak amplitude speed has shifted from 39 000 to 43 000 rpm due to greater 
damping in the squeeze film at the higher eccentricity. The rigidly supported rotor, in 
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contrast, has an amplitude of 50 times the mass eccentricity at the first critical speed. 
Thus the squeeze-film support reduces the rotor amplitude by a factor of 16. 

Figure 17 shows that bearing force at the first critical speed is only 51 newtons 
(11 lb) with the flexible support and that with a rigid support it is 1900 newtons (420 lb), 
or 36 times greater. Similar load reductions occur at the higher critical speeds. At 
28 000 rpm the load for the rigidly supported bearing is 26 000 newtons (5800 lb). This 
far exceeds the load rating of 1530 newtons (344 lb) for the 204 series bearing. With 
the squeeze-film support, however, the load is only 480 newtons (108 lb), well within 
the bearing capacity. 

Figure 18 shows the resultant damper amplitudes for the flexibly supported rotor. 
As predicted for a single -mass rotor, damper amplitudes are approximately equal to 
the mass eccentricity up to twice the first critical speed. At higher speeds, damper 
eccentricities increase. For a large unbalance the increase is about 70 percent; for 
the small unbalance it is nearly a factor of 5 greater. 

Though the damper was sized only for the first critical speed, the results show that 
amplitudes and forces at the higher critical speeds are also reduced substantially. Thus 
single-mass rotor data are useful, not only for multimass rotors operating below the 
second critical speed, but for this rotors also at higher speeds. However, for rotors 
operating through several critical speeds, the authors recommend that a rotor response 
analysis be used after the design of the damper to determine damper and rotor perform- 
ance at the higher critical speeds. 

Rotor response was also calculated with the unbalance concentrated at the center 
mass, rather than distributed over the five masses. This resulted in a much larger 
unbalance loading than with distributed unbalance. Consequently, the squeeze-film 
damper was overloaded; rotor amplitudes and bearing loads were approximately double 
those for rigidly supported bearings. The point to be noted is that an improperly de- 
signed squeeze -film damper (inadequate clearance, etc. ) can be worse than a rigid 
bearing support. Reference 10 also illustrates this phenomenon. 


SUMMARY OF DESIGN PROCEDURE 

(1) For a symmetric multimass rotor, determine an equivalent single mass M 2 
from 


m 2 


*2 


u> 


2 

cr 
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where Kg is the rotor -bearing stiffness at midspan (given by eq. (5)), and u> cr is the 
first bending critical speed for the multimass rotor on rigidly supported bearings. 

(2) Calculate the mass ratio M = Mj/Mg for the total bearing mass Mj to be used 
and the equivalent rotor mass Mg. 

(3) From figure 8 determine the optimum damping ratio B. Determine the absolute 
support damping B required from BjsBX Bg. For lightly damped bearings, the 
value of effective bearing system damping Bg may be estimated from Kg/lOo; . 

(4) Determine the absolute support stiffness K from 

Kj = K X Kg 

For a tuned system, which produces near -minimum rotor amplitude, K = M. Figure 8 
and the design examples of this report assume that a tuned system is to be used. 

(5) Assume a damper eccentricity ratio e that will not make the overall stiffness 
of the support too large, and from figure 14 determine values of dimensionless stiffness 
K d and damping B d _ Generally, the maximum value of e' should be less than 0. 4. 

(6) From the actual rotor mass and the maximum anticipated unbalance, calculate 
the mass eccentricity from 


M 


U 

Mg 


For an optimally damped system (as in fig. 15), the maximum support amplitude is ap- 
proximately equal to the mass eccentricity, that is, aj/e = 1. 

(7) The damper clearance can now be calculated from 



e 


(8) For a circumferentially grooved damper, the damper half length is determined 
from the following: 


L = 



The values of total support damping Bj, the squeeze damping B d , and the clearance C 
all have been previously determined. The damper radius R will usually be dictated by 
the size of the rolling or sliding bearing to be used. The viscosity is that of the 
bearing lubricating oil. 
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(9) Calculate the stiffness of each centering spring from 

K c - 1/2 (k, - 4K d M^ cr RL 3 /c?) 

SUMMARY OF RESULTS 

Theoretical data for a single-mass rotor were used to determine flexible support 
properties (stiffness and damping) to attenuate rotor amplitudes and bearing loads for a 
multimass rotor operating through the first bending critical speed. An equivalent single 
mass for the multimass rotor was calculated from the rotor first critical speed (deter- 
mined from a critical speed computer program) and the rotor shaft stiffness. A 
squeeze-film damper support was then designed to provide the required damping at the 
assumed unbalance conditions. Analytical rotor response results showed that: 

1. The squeeze -film damper successfully attenuated rotor amplitudes and bearing 
loads at the first critical speed. Rotor midspan amplitude was reduced by a factor of 16, 
and bearing load was reduced by a factor of 36 compared with an identical rotor with 
rigidly supported bearings. 

2. Amplitudes and forces at higher critical speeds were also reduced substantially. 

3. With unbalance less than the design value, amplitude and forces were also well 
controlled. However, with unbalances much greater than the design value, amplitudes 
and forces were larger than with rigidly supported bearings. 

4. Bearing loads are well under permissible values for the flexibly supported rotor. 
With rigid supports, bearing forces are very high near rotor critical speeds, resulting 
in drastically shortened bearing lives. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 7, 1974, 

505-04. 
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Figure 2. - Bearing force and rotor amplitude for single-mass flexible rotor in rigidly 
mounted bearings. Amplification factor, 10. 
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Figure 9. - Dimensionless force transmitted to foundation as function of speed ratio for various values of support clamping (from ref. 3). 
Stiffness ratio, 1 ; mass ratio, 1; amplification factor, 10. 
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Figure 10. - Dimensionless force transmitted to bearings as function of speed ratio for various values of support damping (from ref. 31, 
Stiffness ratio, 1 ; mass ratio, 1; amplification factor, 10. 
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Figure 13. - Squeeze-film damper bearing in fixed and rotating coordinate systems. 
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Figure M. - Dimensionless damping (eq. (20)) and 
stiffness <eq, ( 19 )) as functions of eccentricity ratio. 
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